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Abstract 

We study a two-dimensional system of two Coulombically interacting electrons in an external 
harmonic confining potential. More precisely, we present calculations for the singlet ground-state 
of the system. We explain the nature of the degeneracy in the spectrum of the reduced density 
matrix and provide detailed results for the dependencies of the entanglement on the interaction 
strength. Among other features, it is indicated that in the limit of an infinitely strong interaction, 
the angular and radial correlations are asymptotically independent of each other. Appearing in 
this limit, the pure radial correlation effects are quantitatively investigated for the first time. 
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1. Introduction 

Studies of quantum systems trapped in confining potentials have attracted attention because 
of their possible applications in quantum information technology The simplest candidate 
for studying entanglement properties is the system composed of two Coulombically interacting 
electrons in an external harmonic potential, which serves well as a model of a quantum dot (QD)0]. 
In this Letter we consider a two-dimensional (2D) system with a Hamiltonian 

^^2m, 2 e, r2-ri ' ^ ^ 

where m^, is the effective electron mass and e^, is the effective dielectric constant. This system 
was considered in various theoretical contexts [1, d, d, 0, 0| but to the best of our knowledge its 
entanglement properties have not been studied extensively in the literature. 

Since the Hamiltonian ([1]) is spin independent its stationary states possess the form 

^(Ci,C2) =V'(ri,r2)^>(si,S2), (2) 



where g denotes a spin function. After the scaling r i— )• y ;^i7jr, E i— )• the Schrodinger equation 
takes the form 

/?V(ri,r2) =£;V(ri,r2), (3) 
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^ = E[4^r,+2r.^] + -^, (4) 



where the dimensionless couphng g = represents the ratio of Coulomb interaction strength, 

which is related to the so-called Wigner parameter of Ref . 0] by g = ■ 

Introducing the centre-of-mass (cm.) R = ^(ri +r2) and relative (rel.) r = r2 — ri coordinates, 
the Hamiltonian Q is separated into H = H^ + H^ , where the cm. Hamiltonian = — V^/4 + 
AR^ is exactly solvable. The problem is thus reduced to the Schrodinger equation for the rel. 
motion 

described by the Hamiltonian = — + r"^ + ^- For the sake of simplicity we restrict our 
investigation to the study of the singlet ground-state. The total wavefunction is given by 

*(Ci,C2) = ^{r,R)Qs{si,S2) = 
2 _2r2 1 M'"(r) 

where the unknown function (/o°°[''^'"]^^^ = 1) is the lowest-energy eigenfunction of the rel. radial 
equation 

[-|^-^ + r2 + fK(0 = eV(r), (6) 

and Qs = -^(a(l)/?(2) — a(2)/3(l)) is the spin singlet function, where a(/3) denotes the up (down), t 

(I), spin. We note that r = + — 2rir2COs{ip2 — ^i) and R = \/rf + r2 + 2rir2Cos{ip2 — V3i)/2, 
so that in a correlated description ijj depends explicitly on ri, r2 and the inter-electronic angle co- 
ordinate 9 = ip2 — ^pi- 

In this Letter we will analyse the entanglement contained in the state |^') ([5]). Special attention 
is given to performing both the partial-wave expansion and the Schmidt decomposition of tp. We 
make a detailed investigation of the nature of the degeneracy in the spectrum of the reduced 
density matrix. We discuss the effect of interaction on various characteristics such as the collective 
occupancies, the participation ratio, and the Slater rank as well, over a wide range of values of 
g. In particular, we show that only the Schmidt modes with nodeless radial functions contribute 
significantly to the total wavefunction. Importantly, we give an indication that in the 5 — )• oo 
(oj — > 0) limit the global correlation is fully separable into purely angular and radial correlations. 
Moreover, in the regime of — t- 00 we investigate for the first time the role played by the radial 
correlation on the electron pair. 

The Letter is arranged as follows. In Section [2] we discuss the characteristics of the system 
under consideration. Section [3] is devoted to the results and a brief summary of our conclusions is 
provided in Section HI 

2. Entanglement characteristics 

A tool to investigate two-body correlations is the reduced density matrix (RDM) defined as @] 

PrediCX')=tra\^)m, (7) 



where tr denotes a trace taken over one of the electrons, and ^ is the total two-electron wavefunc- 
tion. As ^ is given by Eq. Q, the RDM factors into spatial and spin components 



PrediCC ) = P(r,r )psis,s ), 



(8) 



where the spatial RDM is given by p{r,r') = J ip* {r,r2)ip{T' ,r2)dr2- In particular, for the singlet 
states, the spin part of the RDM is a matrix ps = diag{l/2, 1/2). It is well known that for identical 
fermions the bi-partite pure state I*!*) can be expressed as a combination of the Slater determinants 
made out of one-particle spin-orbitals in which the RDM ([7]) is diagonal The number of non- 
zero expansion coefficients in the Slater decomposition is called the Slater rank (SR) and a pure 
fermion state is entangled if, and only if, its Slater rank is larger than 1. The entanglement depends 
on the whole spectrum of the RDM and many ways of measuring its amount have been developed 
in the literature ld.llll.ll2l.ll3l|. It is interesting to point out that most of the popular entanglement 
measures are functions of the purity of the RDM, trp^^^. The best known are the participation 



ratio 



id ] and the linear entropy ll|. Here we consider the participation ratio. 



redi ' 



(9) 



which approximately counts the number of Slater orbitals actively involved in the Slater expansion 
of ^. The effective Slater rank SR is defined by SR = R/2, since the Slater determinant is 
made up of two different spin orbitals. The larger the value of SR, the higher the entanglement 
or quantum correlation. It is worth stressing at this point that the linear correlation entropy L, 
which is the popular measure of entanglement in pure states 14, 15, [l^, 17, 18 1, is related to R 
via L = 1 - 1/R. 



2.1. The Slater- Schmidt decomposition 

The Slater decomposition of the total wavefunction of ([2]) can easily be inferred from the 
Schmidt decomposition of the spatial part ^. In a previous work |19l |. we have shown that for any 
2D system the real spatial wavefunction that depends only on r and R, il:[r,R) = ij){ri,r2,cos6), 
can be expanded as 



m=— 00...00 



where Arr, = A. 



/rir2 ^/2^^ \/2tt 



(10) 



m — — OO . . . oo 

s — 0,00 



Am{n,r2) = ^^ r r 4,{r,R)e-'"'^^e'^^'d^id^2 
2vr Jq Jo 



2n 



\Jr\r2 \ ip{ri,r2,cos9)cos{m 
'o 



and 



l^s,mXs,m 



(12) 



s=0 
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The radial orbitals Xs,m('') and the coefficients Ks^m are determined by the integral equation 

{r2)dr2 = (n), (14) 





and the family {us^m{T^)} forms a complete and orthogonal set {{usmlu^'^') = /g°° Jq^ ru*^u^i^idrdip 
= ^mm lo" 'X^smXs'm^'^^ ^mm'^ss')^ which means that the latter expansion appearing in Eq. ([lOD 
represents the Schmidt decomposition of ij:. Since the components Am and A^m are the same, 
Eq. ()10p can be written in a partial form 

^{r, R) = MllI^ + y Ai{rur2)cosm ^ (,5^ 

and Xs,m = Xs,-m, i^s,m = i^s.-m- The latter means that the Schmidt decomposition possesses a 
degenerate spectrum, and so it fails to be unique 0]. For the sake of completeness, we give below 
another from Eq. (jlOp form of the Schmidt decomposition of il^. To begin with, we extend the 
results of Ref. [9| to the case with more than one point of double degeneracy. Accordingly, from 
the orbitals Us^i and Us-i (fTHj) (/ > 0), that correspond to the same Schmidt coefficient Ks,h we 
define the new orbitals as 

Us,i{r) + Us,^l{r) .Ms „i(r) - u^j(r) 
^s,lW = ,VsA^) = « ^ — ■ — ' 

be obtained analytically. After some tedious algebra, one arrives at 

j^s,m = — 7= ^,iJsA^) = — 7= (16) 



and for I = as Us o = '^z''" L = Uso, Vgo = 0. The angular parts of the new Schmidt modes can 



The family {i's,u'^s,i} forms a complete and orthonormal set, since one can easily check that 

Wsilt^s'i') ~ ^ss'^ii'^ i'^s,i\fJs'i') ~ ^ss'^ii'^ ('^s.'l^s'z') ~ 0- terms of the new orbitals, the wave- 
function takes the form 

oo 

i'{'r,R) = ^ '^s,iWs,i{^^)^sA^'2) + Vs^i{ri)vs^i{r2)], (17) 

l,s=0 

which yields a Schmidt form different from Eq. ()10p . Consequently, there exist also two different 
forms of the Slater expansion of (0) , which can easily be inferred from the Schmidt decompositions 
(UnD and (HZl). 

2.2. The partial expansion of the RDM 

Now using Eq. (|10p we turn the spatial RDM into the partial form 

p{r,r') = Poir,r^^yPliry)cos{W')^ 



27rv rr j^^^ vrr ^ 



wherein O' = — ip ^ and pi (the /-matrix) is given by 

/>oo 

Pl{r,r')= / Ai{r,r2)Ai{r' ,r2)dr2. (19) 
Jo 

It can easily be derived that the /-matrix has the expansion 

oo 



s=0 



which represents nothing else but its Schmidt decomposition. The eigenequation of the /-matrix is 
thus of the form 

poo ^ ^ ^ 

Pl{r,r)xs,l{^)dr = X^iXs,i{r). (21) 





Obviously, the spatial RDM (jlSp is diagonal in the basis of Schmidt modes (spatial natural orbitals) 
and its eigenvalues A^^; (occupancies) are related to Ks^i via the formula A^^^ = k^^. Moreover, all 
occupancies but the ones with / = are doubly degenerate so that the normalization condition 
gives X^^o A,,o + 2 E£i E^o ^s,i = 1- 

The quantity of interest is the sum of the occupancies corresponding to the /-matrix iji = 
Yl^o^s,i- Throughout this Letter we will refer to {r//} as the collective occupancies. They can 
simply be calculated by r]i = Trpi = pi{r,r)dr = which is crucial since for their 

determination the numerical diagonalizations of the matrices pi are not needed. It is clear that, 
if r]i is vanishingly small compared with the remaining collective occupancies, we may expect that 
the /-partial wave in Eq. (jl5p contributes little. Moreover, 770 might be seen as the fraction of 
electrons with zero angular momentum, and 2r]i (/ > 0), as the fraction of electrons with fixed 
opposite angular momenta {hi and —M). This is particularly easy to understand when rewriting 
the spatial wavefunction as a combination of symmetrized products of the modes u^^i, Ug-i (jl3p q 

2.3. The participation ratio 

The purity of the RDM, trp"^^^, of the singlet state separates in ^trp"^, where ^ is related to the 
spin degree of freedom. Moreover, in the case of (jl8p . trp"^ can be decomposed as 

00 

trp^ = trpl + 2Y,trpl (22) 
1=1 

where 

poo 

trpi= / P^(r,r)dr, (23) 
JO 

with 



00 



Pl{n,r2)= Pi{ri,r3)pi{r3,r2)dr3. (24) 
Jo 

Accordingly, in terms of ()22p the participation ratio of the singlet state can be written as 

1 °° 

B.= [-trpl + Y,trpf]-\ (25) 



1=1 



^ ?/'(r, R) = Zl^o i^s.oUs,o{ri)us.o{r2) + I]^o,i=i V^f^s,i[^ius,i{ri)us,-i{r2) + iisj (r2)us,-i (ri))] 
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From a practical point of view the above expansion appears to be a convenient tool to determine 
the number Ni of partial waves that are needed to capture most of the electrons correlation. It 
seems to be reasonable to define Ni as the minimal number of terms in the sum in (j25p at which 
an approximate value for R has its integer part equal to the integer part of the exact value of R. 

Clearly, the effective Schmidt number of non-zero coefficients in the Schmidt decomposition of 
Ai, Eq. (|12p . is determined by the participation ratio of the normalized Z-matrix, 

\ Pi(r,r') , ^ 

pi{r,r) = ^^^^^, 26 
Trpi 

{Trpi = 1), i.e. [trpi^]-\ 



3. Numerical results 

To start with our analysis, we need the solutions of dH]). As was shown by Taut jj] for a countably 
infinite set of g values the closed-form wavefuntions can be derived, g = \pl is the smallest value 
at which an exact ground-state wavefunction is known. Below this value no analytical solutions to 
the lowest eigenvalue of ([6]) are known and we have determined them numerically by the Rayleigh- 
Ritz method. As was already discussed, a qualitative insight into the character of the partial-wave 
expansion of "0 can be provided by examining the values of the collective occupancies. In Fig. 
[1] the numerically determined behaviour of rn is presented for / = — 5, where to highlight the 
changes with g a logarithmic coordinate has been used. We may notice that at < 1 the collective 




Figure 1: The dependence of the six lowest collective occupancies r]i (/ = — 5) of the lowest singlet state on Ing. 
The inset highlights their changes for large values of g. 

occupancy with / = is about 1. In that situation the angular correlation effects may be neglected 
since one can expect that the terms higher than / = in ()15p contribute very little. When g 
increases, more and more collective occupancies with different I become substantial, which reflects 
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the fact that the larger is the value of g, the larger is the number of partial waves actively involved 
in the sum in ()15p . Interestingly enough, all the collective occupancies except for the largest one 
(770) exhibit a local maximum. The larger is I, the larger is the value of g at which the maximum 
occurs. For large enough values of g the situation changes essentially, namely the clustering of 
collective occupancies is visible. It may have a general interest to note at this point that at 5 < 7.4 
(In 7.4 ~ 2) the collective occupancies with / > 1 are vanishingly small, whereas a Wigner molecule 
is expected to occur when g exceeds a value gcr — 2.8 (See Fig. 5 in Ref. 0], remember that 
g = Rui^/2). Thus, in the incipient Wigner molecule regime, the angular correlations come mainly 
from the first angular level (/ = !). Further, in the range of g in which = rji is substantial, 

we analyse the dependence of the effective Schmidt number for A/, Eq. (|12|) . on g. Our numerical 
results are presented in Fig. [2] in terms of the inverse participation ratio, trpi'^= Qi, for / = — 3 
as a function of lug. The results obtained indicate nothing but that those partial waves that 
contribute considerably to ip have their radial parts Ai being almost uncorrelated functions. That 
is, Ai{ri,r2) ~ Xo,i{'''i)Xo,ii^2), as we have verified by solving Eq. (fT^ through a discretization 
technique. In other words, we have the result that only the Schmidt orbitals with (s = 0, /) in ()17p 
are important and, in consequence, the total wavefunction approaches the form 

*(Ci,C2)^^(a(l)/3(2)-a(2)/3(l))x 

Ni-l 

XI '^o,iWo,i{^i)i^o,i{^2) + vo,i{ri)vo^i{r2)], (27) 
1=0 

where Ni is the number of substantial partial waves. The inspection of (j27p yields that it constitutes 
a sum of 2A''; — 1 Slater determinants (here we recall that t;o,o = 0). In the limit of g — )• 00, the 
values of {^i}i=o saturate to the same value ^1^°°^ ~ 0.9634, which reflects the fact that at (7 —t- 00 
the occupancies are asymptotically independent of the angular momentum i.e., A^^°° = A^^°° = 

... = Ai°°'' (/Oq~^°° = pf^°° = ... = p''°°^). We find this consistent with the properties of the 3D 
Hooke's atom that were discussed in [2^. What follows further from the above is that in the limit 
g ^ 00, the radial components {Ai}i=q can differ from each other only by their signs. We have 
verified that Af^°° = {—lY A^°°\ Substitution of this into Eq. (jlSp gives 

(r,i?) = —^^i^ + 'Ei-iycosiW)]. (28) 

Because (— l)'cos(/0) = cos[l{6 — vr)], the suitable part of Eq. (pHj) can be recognized as a Dirac 
delta distribution centred at vr i.e. The angular correlations carry thus the electrons 

strictly at opposite sides of the centre of the trap as (7 — )• 00, which corresponds to a perfect linear 
Wigner molecule. In this limit the potential energy dominates over the electron kinetic energy and 
the electrons localize to the classical positions = |(f )^) ^ = '^i but always there exist quantum 
fluctuations with respect to their radial coordinates. The electrons oscillate about the equilibrium 
geometry in accordance with a harmonic approximation model 0, [E; S]- At g — )• 00 there is no 
fluctuation in 6, but the angular correlations are extremely strong due to an arbitrary orientation 
of the molecular axis, which coincides with the fact that the lowest-energy classical configuration is 



■^A Fourier series expansion of 6(x — s), gives 5{x — s) = + ^ cos[x{l — s)][2^ 
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infinitely degenerate with respect to rotations around the center of the trap. When it comes to pure 
radial correlations appearing in the asymptotic regime they are expected to be very poor, due to the 
closeness of the value of 0^°°) to unity (0^°°) ~ 0.9634). More precisely, we have found numerically 
that the ratio of X^^/X^f°^ is about 54, 2883, 154804 for s = 1,2,3, respectively, and it tends 
to rapidly increase with further increasing s. In other words, the radial component appearing in 
Eq. ()28p can indeed be considered as almost a product function, which in turn implies low pure 
radial correlation effects. We stress that the limiting values of the considered quantities have been 
obtained with the use of the wavefunction in the harmonic approximation, which is valid as (7 — )• oo. 
Before going further we want to stress that in none of the well known papers concerning the system 
under consideration has Eq. (|28|) appeared. To gain further understanding we next quantify the 
degree of entanglement by calculating the effective Slater rank SR = R/2. An approximate value 
for R can be obtained with the help of Eq. (j25p by successively increasing the number of terms 
in the sum appearing in it until the result becomes stable to a desired accuracy. In Fig. [3] the 
behaviour of SR determined by the procedure described above is shown. We have checked that 

1.00 
0.99 
0.98 
0.97 



0.96 



2 4 6 8 

Ing 




Figure 2: The behaviour of ili for I — — 3 as a, function of Ing. The first vertical line (on the left) corresponds to 
the value of g below which the partial wave with I = 1 hardly contributes to -0 at all. The second line concerns the 
same but for I = 2, and so on. The locations of the lines were roughly inferred from Fig. [T] 

the results for the linear entropy calculated with Eq. ()25p as L = 1 — 1 /R completely agree with 
the results presented in Ref. [2l|], where they were calculated in a different way. The circles in 
Fig. [3] represent the values of the Slater rank obtained from our earlier finding SR= 2A''; — 1 with 
Ni being determined as was discussed at the end of subsection 12.31 As one can see, the behaviour 
of SR coincides fairy well with that of SR, which justifies the approximation (|27p . 



4. Summary 



In conclusion, we have investigated the nature of the degeneracy in the spectrum of the spatial 
RDM of the two-particle wavefunction that is a function of distances r and R only. We have devel- 
oped the results of Ref. 19, 2^ by providing the Schmidt decomposition in the single real particle 




2 4 6 8 10 



Ing 

Figure 3: The effective Slater rank SR = R/2 as a function of Ing. The circles mark the results obtained by 
SR= 27V, - 1. 

basis in parabolic coordinates. In addition, we carried out a comprehensive study of the singlet 
ground-state of two interacting electrons confined in an isotropic harmonic potential. As a general 
trend we found that as g increases, the number of partial waves that are needed to capture most of 
the electrons' correlation increases as well. The Slater rank grows monotonically with increasing 
g and the bulk of the entanglement is mainly manifested in the angular variables. At (7 — )• oo the 
spatial wavefunction factorizes asymptotically into a product of radial and angular components, 
of which the former is individually little correlated, but the latter is maximally correlated. Our 
calculations have shown that only at 5 < 1 can the total wavefunction be well approximated by 
a single Slater determinant. In other words, in this regime the state under consideration can be 
regarded as weakly entangled. 

It would be interesting to apply the tools presented in this Letter to explore the entanglement 
properties of two-particle systems with other interactions between the particles. This will be the 
topic of our further studies. 
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